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SUMMARY 


This  technical  note  outlines  the  MICE/MEGS/MELT  -  SCENARIO 
Plasma  Outputs  Interface  specification,  Version  1.031. 

High-altitude  nuclear  phenomenology  codes  at  the  first- 
principles  level  (MICE/MEGS/MELT  codes)  and  at  the  engineering-level 
(SCENARIO  multiburst  code)  generate  computer  files  which  contain  a 
variety  of  descriptors  of  the  disturbed  environment  resulting  from  high- 
altitude  nuclear  detonations.  Efficient  and  interchangeable  use  of 
results  from  thes*  codes  requires  a  common  interface  specification  with 
sufficient  generality  to  permit  an  end-user  to  accommodate  basic 
differences  between  the  codes  and  their  outputs.  This  document 
represents  the  latest  version  of  the  standard  file  format.  Publication 
coincides  with  distribution  of  SCENARIO,  Version  4.0,  which  utilizes  this 
Interface. 
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SECTION  1 


INTRODUCTION 


This  document  describes  a  revised  Standard  Plasma  Outputs 
Interface  format  which  is  common  to  the  SCENARIO,  MICE/MEGS/MELT  family 
of  nuclear  phenomenology  codes  sponsored  by  the  Defense  Nuclear  Agency. 
The  purpose  of  a  standardized  format  is  to  enable  easy  transfer  of 
information  between  nuclear  phenomenology  codes  and  post-processor  / 
end-user  codes. 

In  this  context,  post-processor  codes  augment  basic 
environmental  descriptors  produced  by  the  original  phenomenology 
simulation,  be  it  from  SCENARIO  or  MICE/MELT.  End-user  codes  typically 
use  basic  and/or  post-processor-augmented  phenomenology  outputs  to 
compute  system  performance  parameters  for  operation  in  a  nuclear 
environment. 

For  example,  the  PRPSIM  post-processor  code1  directly  uses 
nuclear  phenomenology  code  outputs  (e.g.,  SCENARIO  or  MICE/MELT  results) 
to  provide  performance  parameters  associated  with  radio  frequency 
communication  links.  An  alternative  operating  mode  directs  basic 
phenomenology  outputs  (in  Standard  Format)  to  an  intermediate 
post-processor.  The  intermediate  post-processor  i)  calculates  additional 
phenomenological  quantities  (e.g.,  structure  descriptors  such  as  outer 
scale,  freezing  scale,  and  spectral  indices  associated  with  a  power 
spectral  density)  and  ii)  writes  in  Standard  Format  the  original 
description  of  the  burst-disturbed  environment  plus  its  additional 
quantities.  Because  the  Standard  Format  is  preserved  in  the  process, 
PRPSIM  and  other  post-processor  /  end-user  codes  can  read  and  process  the 
augmented  file  as  well. 

Although  basic  objectives  of  the  Standard  Plasma  Outputs 
Interface  had  been  accomplished  by  earlier  versions  of  the  format,  recent 
systems  requirements  for  additional  phenomenology  information,  coupled 
with  improved  capabilities  for  predicting  underlying  nuclear  burst 
effects,  have  forced  revision. 

Older  versions  of  the  format  made  allowances  for  future 
extension  to  include  information  about  infrared  environments  from  nuclear 
bursts.  However,  the  infrared  problem  has  proven  sufficiently  complex 
that  it  is  desirable  to  separate  MHD/plasma  descriptors  (electron 
density,  ion  species  densities,  neutral  species  densities,  velocities, 
etc.)  from  information  about  infrared  environments.  The  word  "plasma" 
has  been  inserted  into  the  title  of  this  interface  specification  to 
emphasize  the  separation.  A  parallel  document  describes  the  standardized 
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format  for  infrared  information.2 

In  many  areas,  organization  of  the  Interface  has  been  preserved 
without  change.  To  accommodate  new  approaches  to  plasma  structure 
modeling  and  extended  descriptions  of  the  power-spectral -density  of 
plasma  irregularities,  a  variety  of  changes  to  the  record  structure  of 
the  Standard  Plasma  Outputs  Interface  have  been  made.  These 
rearrangements  were  motivated  by  DMA  requirements  plus  a  desire  to 
minimize  future  changes  as  structure  specifications  become  more  refined. 


To  integrate  new  structure  information  into  the  Interface, 
structure  parameters  have  been  moved  to  a  new  location  within  the 
Environment  Tube  Records.  (Note:  terminology  is  summarized  in  the 
Glossary.)  Environment  Tube  Records  hold  environmental  descriptors 
(electron  density,  velocity,  temperature,  etc.)  on  a  cell  by  cell  basis. 
The  new  arrangement  situates  structure  descriptors  for  ionized  material 
in  a  new  set  of  (optional)  cell  quantities  located  immediately  after  the 
"usual"  quantities  in  the  Environment  Tube  Records.  Structure  quantities 
are  optional  since  not  all  high-altitude  simulations  will  generate 
structure  information.  The  presence  or  absence  of  these  quantities  is 
indicated  by  the  value  of  the  ISTRIR  parameter  in  the  Header  Record.  The 
number  of  structure  quantities  per  cell,  NSTRQ,  is  also  a  new  parameter 
in  the  Header  Record. 


Another  change  of  note  involves  addition  of  neutral  fluid 
velocity  information  to  the  Environment  Tube  Records.  This  was  done  to 
provide  a  basis  for  post-processing  with  structure  modeling  techniques 
which  use  the  differential  velocity  (slip  velocity)  between  the 
ion-electron  plasma  and  the  neutral  fluid. 

Users  of  this  Interface  should  pay  particular  attention  to 
ordering  and  content  of  the  revised  Environment  Tube  Records.  Revisions 
discussed  above  have  forced  a  number  of  significant  changes.  Each  has 
been  flagged  by  an  asterisk  (*)  in  the  margin.  A  few  typographical 
errors  from  earlier  versions  of  the  Outputs  Interface  specification  have 
also  been  repaired. 
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SECTION  2 


INTERFACE  SPECIFICATION 


Numerical  outputs  from  MICE/MELT  or  SCENARIO  phenomenology 
simulations  are  generated  at  specific  simulation  times  selected  by  the 
operator  of  the  code.  Information  written  to  the  program’s  output  file 
is  grouped  into  discrete  clusters  of  records.  Each  cluster  represents 
environmental  conditions  at  a  particular  phase  in  the  simulation.  For 
purposes  of  discussion  below,  each  cluster  of  records  will  be  termed  a 
code  output  block.  In  the  phenomenology  community,  code  output  blocks 
are  sometimes  called  "dumps"  or  "write-outs". 

It  is  important  to  recognize  that  within  an  output  file  for  a 
simulation,  any  code  output  block  may  differ  in  important  ways  from  other 
code  output  blocks  (i.e.,  other  simulation  times).  For  example,  changes 
to  computational  grids  might  occur  during  a  simulation,  altering  numbers, 
locations,  and  dimensions  of  computational  cells.  Changes  of  this  sort 
could  be  triggered  by  a  need  to  alter  grid  resolution  as  the  simulation 
evolves.  Numerous  other  shifts  of  basic  simulation  parameters  are 
possible.  The  message:  one  cannot  assume  basic  simulation  parameters 
remain  constant.  The  Header  Record  for  each  code  output  block  should 
always  be  read. 

The  sub-sections  which  follow  detail  the  structure  of  an 
individual  block  of  code  output.  The  sub-sections  are  organized  to 
reflect  the  sequential  record  structure  found  in  actual  code  output. 
FORTRAN  syntax  is  used  for  mathematical  expressions.  FORTRAN  naming 
conventions  for  distinguishing  Real  variables  from  Integers  are  followed. 
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2.1  HEADER  RECORD. 


Variable 

Description 

IDENT 

Alphanumeric  identification  of  simulation. 

(80  characters) 

SOURCE 

Alphanumeric  identification  of  computer  code 
used  to  generate  simulation.  (80  characters) 

CONTACT 

Alphanumeric  identifier  of  appropriate  people 
to  contact  for  resolution  of  phenomenology 
questions.  (120  characters) 

COMMENT 

General  alphanumeric  comments  field. 

(80  characters) 

RUNDATE 

Alphanumeric  field  containing  time  and  date 
of  simulation.  (20  characters) 

(*) 

ISOURCE 

Integer  simulation  source  code  ID: 

1  =  MICE 

2  -  MEGS 

3  =  MELT 

4  =  SCENARIO 

VERSION 

Real  variable  denoting  version  number  of 

Plasma  Outputs  Interface  specification. 

Current  version  is  1.031,  implemented  6/89. 

(*) 

TSIM 

Simulation  time  (seconds)  for  this  block  of 
code  output. 

ITIME 

Simulation  clock  flag: 

0  =  Simulation  time  relative  to  an  arbitrary 
reference  (usually  the  first  burst). 

1  =  Simulation  time  is  GMT. 


NB  Number  of  bursts  detonated  prior  to  the  writing 

of  this  block  of  code  output. 

NRECS  Number  of  records  (including  this  one)  in 

the  current  block  of  code  output. 


4 


IGEOMR 


ISYMR 


IFGRID 


NATUBES 


Geometry  flag: 

1  =  MICE  3-D  modified  Geocentric  Spherical 

coordinates . 

2  =  Geomagnetic  Dipole  coordinates 

with  aligned  a  array. 

3  =  Geomagnetic  Dipole  coordinates 

with  staggered  a  array. 

Grid  symmetry  flag: 

0  =  No  symmetry 

1  =  Magnetic  east-west  reflection 

symmetry  about  <f>=0  (or  Y=0  in 
MICE  output). 

2  =  North-south  reflection  symmetry 

about  magnetic  equator  (or  X=0 
in  MICE  output) . 

3  =  Both  north-south  and  east-west 

reflection  symmetries. 

4  =  Special  symmetry  —  consult  MRC 

phenomenologists. 

Full  grid  output  flag: 

0  =  Output  block  represents  full 

computational  grid.  All  environment 
tubes  -  both  active  and  ambient. 

Total  number  of  tubes,  active  plus 
ambient,  is  JMXB*JMXP.  Number  of  active 
tubes  is  NATUBES. 

1  =  Output  block  contains  active  (non¬ 
ambient)  environment  tubes  only.  There 
are  NATUBES  active  environment  tubes. 
Ordering  of  active  tubes  within  the 
grid  is  as  specified  by  the  Active/ 
Inactive  Environment  Tube  Record. 
Criteria  used  to  determine  which  tubes 
are  active  depend  on  the  particular 
code  (ISOURCE). 

Number  of  active  environment  tubes  in  this 
code  output  block. 
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NBUFF 


JMXA 

JMXB 

JMXP 

NCELLS 

NQPCEL 

NSTRQ 

LMXI 

LMXN 

NAMB 


Buffer  size  required  to  read  Environment  Tube 
Records.  Nominal  implementation  is  one 
environment  tube  per  record,  so  minimum  buffer 
size  is 

(JMXA  *  NQPCEL  +  3).  (*) 

To  minimize  interrecord  gaps  on  magnetic  tape, 
buffer  size  may  optionally  be  increased  to 
greater  than  or  equal  to  an  integer  multiple 
of  the  minimum  size.  Consult  MRC  phenomenolo- 
gists  to  arrange  for  larger  buffer  sizes. 

Number  of  grid  cells  in  "vertical"  direction. 
"Vertical"  is  the  Z  direction  in  MICE,  but  is 
the  a  direction  in  MEGS,  MELT,  and  SCENARIO. 

Number  of  grid  cells  in  north-south  direction. 
"North-south"  is  the  X  direction  in  MICE,  but 
is  the  p  direction  in  MEGS,  MELT,  and  SCENARIO. 

Number  of  grid  cells  in  east-west  direction. 
"East-west"  is  the  Y  direction  in  MICE,  but  is 
the  4>  direction  in  MEGS,  MELT,  and  SCENARIO. 

Total  number  of  grid  cells  in  this  code  output 
block. 

If  IFGRID  «  0,  then  NCELLS  =  JMXA*JMXB*JMXP. 

If  IFGRID  =  1,  then  NCELLS  =  NATUBES*JMXA. 

Number  of  quantities  per  cell  in  this  code 
output  block:  (*) 

30  +  NSTRQ  for  MEGS/MICE/MELT; 

18  +  NSTRQ  for  SCENARIO. 

Number  of  striation  quantities  per  cell.  (*) 

Default  value  is  0  (if  ISTRIR  =  0) 
or  5  (if  ISTRIR  >  0). 

Number  of  ion  species  per  cell. 

(Does  not  include  weapon  debris.) 

Number  of  neutral  species  per  cell. 

Number  of  altitudes  at  which  ambient 
atmosphere/ionosphere  quantities  are  specified. 
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ICHEMR 

IDEBRIS 

FBETA 

FPHOT 

IBETAR 

ISTRIR 


Integer  chemistry  model  flag: 

1  =  Standard  SCENARIO  two  reaction 

chemistry  model. 

2  =  Standard  MICE/MELT  chemistry  model. 

3  =  Standard  MICE/MELT  chemistry  model 

plus  D-region  lumped-parameter 
chemistry. 

Weapon  debris  flag: 

0  =  No  debris  density  specified. 

1  =  Weapon  debris  density  specified. 

Fraction  of  debris  mass  density  that  is  a 
beta  emitter. 

Fraction  of  debris  mass  density  that  is  an 
IR/visible  scatterer. 

Beta  flux  record  flag: 

0  =  No  Beta  Flux  Record  written  in  this 
code  output  block. 

1  *  Beta  Flux  Record  written. 

Striation  model  flag:  (*) 

If  ISOURCE  -  1  (MICE)  - 

0  =  No  striation  parameters  initialized, 
no  striation  quantities  written. 

1  =  Local,  linear  striation  analysis; 
striation  quantities  written. 

If  ISOURCE  »  2  (MEGS)  -- 

0  =  No  striation  parameters  initialized, 
no  striation  quantities  written. 

1  =  Local,  linear  striation  analysis; 
striation  quantities  written. 

If  ISOURCE  =  3  (MELT)  -- 

0  »  No  striation  parameters  initialized, 
no  striation  quantities  written. 

1  =  Field-line-averaged  striation  growth 

plus  regrid-variance  accumulation; 
striation  quantities  written. 

2  =  Field-line-averaged  striation  growth 

plus  regrid-variance  accumulation 
plus  striation  convection; 
striation  quantities  written. 

3  -  Non-equipotential  striation  model 

parameters.  (Not  yet  implemented) 
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CLOR 

ELOR 

ISPARE1 

ISPARE2 

ISPARE3 

ISPARE9 


If  ISOURCE  =  4  (SCENARIO)  -- 

0  =  no  striation  parameters  initialized, 
no  striation  quantities  written. 

1  =  standard  SCENARIO  striation 

convection  parameters; 
striation  quantities  written. 

2  =  non-equipotential  striation  model 

parameters;  striation  quantities 
written.  (Not  yet  implemented) 

Magnetic  colatitude  (radians)  of  origin  of 
geocentric  spherical  grid  (MICE  outputs 
only,  zero  otherwise). 

East  magnetic  longitude  (radians)  of 
origin  of  geocentric  spherical  grid 
(MICE  outputs  only,  zero  otherwise). 
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2.2  BURST  HISTORY  RECORD. 


(TIMEB(I),  ZB( I) ,  XB( I) ,  YB(I),  TYB(I) ,  1=1 ,NBA)  (*) 

where: 

TIMEB(I)  =  Detonation  time  of  the  Ith  burst.  (sec) 

(Bursts  listed  chronologically.) 

ZB(I)  =  "Vertical"  coordinate  of  Ith  burst  point. 

Z  coordinate  (altitude)  in  geocentric 

spherical  system  (MICE);  (cm) 

a  coordinate  in  geocentric  dipole 

system  (MEGS,  MELT,  SCENARIO).  (rad) 

XB{ I )  =  "North-south"  coordinate  of  Ith  burst  point. 

X  coordinate  (positive  to  north)  in  geo¬ 
centric  spherical  coordinate  system  (MICE);  (cm) 

f)  coordinate  (positive  toward  magnetic 
equator)  in  geocentric  dipole  system 
(MEGS,  MELT,  SCENARIO).  (rad) 

YB( I)  =  "East-west"  coordinate  of  Ith  burst  point. 

Y  coordinate  (positive  to  west)  in  geo¬ 
centric  spherical  coordinate  system  (MICE);  (cm) 

4>  coordinate  (positive  to  east)  in  geo¬ 
centric  dipole  system  (MEGS,  MELT,  SCENARIO).  (rad) 

TYB(I)  =  Total  Yield  of  Ith  burst.  (MT)(*) 


NBA  =  MAX(NB, 1) 

(See  Header  Record  for  description  of  NB.) 


NOTE:  If  no  detonations  have  occurred  prior  to  the  first  code  output 
block,  then 

NBA  =  1,  TIMEB  *  -1.E30, 

ZB ( 1 )  =  XB ( 1 )  =  YB( 1 )  =  0., 

TYB(l)  =  0. 
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2.3  GRID  GEOMETRY  RECORD. 


(  (  ZKC(I,J),  I-l.JMXA  ),  J-l.JMXB  ), 

(  (DZKC(I.J),  1=1 , JMXA  ),  J-l.JMXB  ), 

(  XKC(I),  1=1,  JMXB  ), 

(DXKC(I),  1=1,  JMXB  ), 

(  YKC(I),  1=1,  JMXP  ), 

(DYKC(I),  1=1,  JMXP  ) 

where: 

ZKC  =  Array  of  cell -centered  "vertical"  coordinates. 

Z  coordinate  (altitude)  in  geocentric  spherical 
coordinate  system  (MICE);  (cm) 

a  coordinate  in  geocentric  magnetic  dipole  coordinate 
system  (MEGS,  MELT,  SCENARIO).  The  a  unit  vector 
is  anti -paral lei  to  the  dipole  geomagnetic  field.  (rad) 

NOTE:  ZKC  is  doubly  subscripted  to  accommodate 
staggered  a  grids. 

DZKC  =  Array  of  cell  widths  associated  with  ZKC.  (cm  or  rad) 

NOTE:  DZKC  is  doubly  subscripted  to 
accommodate  staggered  a  grids. 

NOTE:  For  MICE  only,  altitudes  of  both  the  (*) 

cell  center  and  lower  boundary  of  the  bottom 
cell  are  defined  to  be  the  same.  This 
oddity  results  from  use  of  a  reflective 
bottom  boundary  condition  in  the  MICE  finite 
difference  scheme. 

XKC  =  Array  of  cell -centered  "north-south"  coordinates. 

X  coordinate  (positive  toward  north)  in  geocentric 


spherical  coordinate  system  (MICE);  (cm) 

p  coordinate  (positive  toward  magnetic  equator) 
in  geocentric  dipole  coordinate  system  (MEGS, 

MELT,  SCENARIO).  (rad) 

OXKC  =  Array  of  cell  widths  associated  with  XKC.  (cm  or  rad) 

YKC  =  Array  of  cel  1 -centered  "east-west"  coordinates. 

Y  coordinate  (positive  toward  west)  in  geocentric 
spherical  coordinate  system  (MICE);  (cm) 

4>  coordinate  (positive  toward  east)  in  geocentric 
dipole  coordinate  system  (MEGS,  MELT,  SCENARIO).  (rad) 

OYKC  =  Array  of  cell  widths  associated  with  YKC.  (cm  or  rad) 
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2.4  AMBIENT  ATMOSPHERE/ IONOSPHERE  RECORD. 


If  ISOURCE  *  1,  2,  or  3  (MICE,  MEGS,  or  MELT): 

(Z(I),  ENE(I),  RHON(I),  RHO 1(1),  TEMP(I),  (SPEC(J,I),  J-1,8), 

I-l.NAMB)  (*) 

If  ISOURCE  -  4  (SCENARIO): 

(Z(I),  ENEDAY(I) ,  ENENITE(I),  TEMP(I),  RHON(I),  'I-l.NAMB) 
where: 

Z(I)  =  Altitude  of  Ith  entry  in  ambient  table.  (cm) 

ENE(I)  -  Ambient  electron  density  at  Ith  altitude.  (cm'3) 

ENEDAY(I)  =  Daytime  ambient  electron  density  at  Ith  (cm'3) 

altitude.  (SCENARIO  only.) 

ENENITE(I)  =  Nighttime  ambient  electron  density  at  Ith  (cm*3) 

altitude.  (SCENARIO  only.) 

RHON(I)  *  Ambient  neutral  mass  density  at  Ith  altitude,  (gm/cm3) 

RHOI(I)  -  Ambient  ion  mass  density  at  Ith  altitude.  (gm/cm3) 

TEMP(I)  »  Ambient  temperature  at  Ith  altitude.  (deg  K) 

SPEC (1,1)  =  Ambient  N2  density  at  Ith  altitude.  (cm"3) 

SPEC (2,1)  =  Ambient  02  density  at  Ith  altitude.  (cm'3) 

SPEC (3,1)  -  Ambient  0  density  at  Ith  altitude.  (cm'3) 

SPEC (4 , 1 )  -  Ambient  NO  density  at  Ith  altitude.  (cm'3) 

SPEC (5,1)  =  Ambient  H  density  at  Ith  altitude.  (cm'3) 

SPEC (6,1)  =*  Ambient  0+  density  at  Ith  altitude.  (cm'3) 

SPEC ( 7 , 1 )  -  Ambient  NO"  density  at  Ith  altitude.  (cm'3) 

SPEC (8 , 1 )  -  Ambient  H+  density  at  Ith  altitude.  (cm'3) 

NOTE:  See  Header  Record  for  description  of  NAMB. 
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2.5  ACTIVE/INACTIVE  ENVIRONMENT  TUBE  RECORD. 


( ( I MA P ( I , J) ,  I-l.JMXB),  J-l.JMXP) 
where: 


IMAP(I,J)  *  Index  map  for  active  and  inactive  environment  tubes. 
The  (I,J)  index  corresponds  to  the  (I,J)th  {p,<t>) 
grid  location  where  1  <  I  <  JMXB  and  1  <  J  <  JMXP. 

0  =  (I,J)th  environment  tube  is  inactive. 

IFGRID  parameter  in  Header  Record  controls 
whether  or  not  (I,J)th  environment  tube 
is  written  to  output  block. 

n  *  Positive  integer  which  indexes  active 
environment  tubes.  Active  tubes  are 
written  and  numbered  in  sequential  order 
as  grid  is  scanned  most  rapidly  along 
I  index  (0  direction),  least  rapidly  along 
J  index  (<f>  direction). 
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2.6  ENVIRONMENT  TUBE  RECORDS. 

( I B  ETA , J  PH I , KTUBE ,(( CQBUFF ( NQ , NC , NT ) , 

NQ=1,NQPCEL),  NC=1 , JMXA) ,  NT=1, NTBUFF)  (*) 


where: 

IBETA  =  p  index  of  environment  tube. 

JPHI  =  <t>  index  of  environment  tube. 

KTUBE  =  Active  environment  tube  index. 

0  =  Inactive  environment  tube  (if  written), 
n  =  Active  environment  tube  index  as  in 
IMAP( IBETA, JPHI) . 

CQBUFF  =  Array  used  to  read  cell  quantities  in  environment 
tubes.  Each  tube  record  holds  cell  quantities  for 
one  or  more  complete  environment  tubes.  No  partial 
tubes  are  written.  The  number  of  complete  tubes  in 
a  full  buffer  is  the  integer 

NBUFF/(NQPCEL*JMXA+3)  (*) 

The  first  CQBUFF  index  runs  over  cell  quantities: 

1  <  NQ  s  NQPCEL 

The  second  CQBUFF  index  runs  over  cells  in  one 
environment  tube: 

1  <  NC  <  JMXA 

The  third  CQBUFF  index  runs  over  environment  tubes: 

1  <  NT  <  NTBUFF 

NOTE:  For  NTBUFF  >  1,  the  last  Environment  Tube 

Record  will  generally  contain  fewer  than  NTBUFF 
tubes.  This  occurs  because  the  total  number  of 
tubes  may  not  be  an  integer  multiple  of  NTBUFF. 


CQBUFF ( 1 ,NC,NT)  =  Cell  index  (real  variable)  which  runs  over 

all  cells  (active  plus  inactive)  in  grid. 

This  index  is 

NC  +  JMXA  *  (  IBETA- 1  +  JMXB  *  (  JPHI-1  )) 
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(cm'3)  (*) 


CQBUFF ( 2 , NC , NT ) 


CQBUFF(3,NC,NT) 


CQBUFF ( 4, NC, NT) 


CQBUFF(5,NC,NT) 

CQBUFF(6,NC,NT) 

CQBUFF(7,NC,NT) 


CQBUFF(8,NC,NT) 

CQBUFF(9,NC,NT) 
CQBUFF ( 1 0 , NC , NT ) 


Electron  number  density. 

This  quantity  will  be  the  mean  electron 
density  if  striation  convection  is 
operative;  i.e., 

if  ISOURCE  -  3  and  ISTRIR  =  2 

or  ISOURCE  =  4  and  ISTRIR  *  1 

Partial  time  derivative  of  electron  (cm'3-sec'1) (*) 
density.  This  quantity  will  be  a  mean 
value  if  striation  convection  is  operative. 

(This  quantity  will  be  zero  unless 
specifically  requested  by  the  end  user. 

Not  available  from  standard  SCENARIO  simula¬ 
tions  at  present.) 

Debris  mass  density.  (cm'3)(*) 

Ionized  weapon  debris.  This  quantity  will 

be  a  mean  value  if  striation  convection 

is  operative.  (This  quantity  is  non-zero  only 

if  debris  was  initialized  and  transported  in 

the  simulation.  See  also  IDEBRIS  description 

in  Header  Record.) 

Ion  mass  density.  (gm/cm3) 

Ion  velocity  in  the  "vertical”  direction,  (cm/sec) 

Z  direction  in  geocentric  spherical  grid 
(MICE); 

a  direction  in  geocentric  dipole  grid 
(MEGS,  MELT,  SCENARIO). 

Ion  velocity  in  the  "north-south" 

direction.  (cm/sec) 

X  direction  in  geocentric  spherical  grid 
(MICE); 

p  direction  in  geocentric  dipole  grid 
(MEGS,  MELT,  SCENARIO). 

Ion  velocity  in  the  "east-west”  direction,  (cm/sec) 

Y  direction  in  geocentric  spherical  grid 
(MICE); 

4  direction  in  geocentric  dipole  grid 
(MEGS,  MELT,  SCENARIO). 

Neutral  mass  density.  (gm/cm3) 

Molecular  mass  fraction. 

Neutral  molecule  mass  density  divided  by 
total  neutral  mass  density. 
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(cm/sec) 


CQBUFF(11 ,NC, NT)  =  Neutral  velocity  component  in  the 

"vertical"  direction.  Z  direction  in  a 
geocentric  spherical  grid  (MICE,  MELT, 

MEGS);  a  direction  in  a  geomagnetic 
dipole  grid  SCENARIO). 

CQBUFF ( 1 2 , NC , NT)  =  Neutral  velocity  component  in  the  "north-  (cm/sec) 

south"  direction.  X  direction  in  a  geo¬ 
centric  spherical  grid  (MICE,  MELT,  MEGS); 

$  direction  in  geomagnetic  dipole  grid 
(SCENARIO). 

CQBUFF(13,NC,NT)  =*  Neutral  velocity  component  in  the  "east-  (cm/sec) 

west"  direction.  Y  direction  in  a  geo¬ 
centric  spherical  grid  (MICE,  MELT,  MEGS); 

<t>  direction  in  geomagnetic  dipole  grid 
(SCENARIO). 

CQBUFF(14,NC,NT)  «  "Vertical"  component  of  magnetic  field.  (gauss) 

Z  direction  in  geocentric  spherical 
grid  (MICE); 

a  direction  in  geocentric  dipole 
grid  (MEGS,  MELT,  SCENARIO). 

CQBUFF(I5,NC,NT)  =  "North-south"  component  of  magnetic  field,  (gauss) 

X  direction  in  geocentric  spherical  grid 
(MICE); 

p  direction  in  geocentric  dipole  grid 
(MEGS). 

Zero  for  MELT  and  SCENARIO  by  assumption  of 
ambient  geomagnetic  field. 

CQBUFF(16,NC, NT)  *  "East-west"  component  of  magnetic  field.  (gauss) 

Y  direction  in  geocentric  spherical  grid 
(MICE); 

4>  direction  in  geocentric  dipole  grid 
(MEGS). 

Zero  for  MELT  and  SCENARIO  by  assumption  of 
ambient  geomagnetic  field. 

CQBUFF(17,NC,NT)  =  Electron  temperature.  (deg  K) 

(Plasma  temperature  if  IS0URCE=4 
(SCENARIO).) 

CQBUFF(18,NC,NT)  *  Neutral  temperature.  (deg  K) 


END  of  Standard  SCENARIO  Quantities 
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Extended  Phenomenology  Quantities 


-  The  following  quantities  will  be  present  if  ISOURCE  =  1,  - 

-  2,  or  3  (i.e.,  MICE,  MEGS,  or  MELT,  respectively).  - 

-  If  present,  use  NS  -  18  below.  - 

-  If  ISOURCE  *  4  (SCENARIO),  these  temperature  and  species  - 

-  quantities  are  not  present.  Skip  to  structure  quantities.  - 


CQBUFF(NS+1,NC,NT)  =  N2  vibrational  temperature.  (deg  K) 

CQBUFF(NS+2,NC,NT)  =  Ion  temperature.  (deg  K) 

CQBUFF(NS+3,NC,NT)  =  N+  density.  (cm'3) 

This  quantity  will  be  a  mean  value 
if  striation  convection  is  operative. 

(See  CQBUFF(2,NC,NT) .) 

CQBUFF(NS+4,NC,NT)  =  0*  density.  (cm'3) 

This  quantity  will  be  a  mean  value 
if  striation  convection  is  operative. 

(See  CQBUFF(2,NC,NT) . ) 

CQBUFF (NS+5 . NC , NT)  =  N0+  density.  (cm'3) 

This  quantity  will  be  a  mean  value 
if  striation  convection  is  operative. 

(See  CQBUFF(2,NC,NT) . ) 

CQBUFF (NS+6,NC, NT)  =  H+  density.  (cm'3) 

This  quantity  will  be  a  mean  value 
if  striation  convection  is  operative. 

(See  CQBUFF(2 ,NC,NT) . ) 

CQBUFF (NS+7.NC, NT)  =  N2  density.  (cm'3) 

CQBUFF (NS+8,NC, NT)  =  02  density.  (cm'3) 

CQBUFF (NS+9,NC,NT)  =  NO  density.  (cm'3) 

CQBUFF (NS+10,NC, NT)  =  N  density.  (cm'3) 

CQBUFF (NS+ 11 ,NC,NT)  =  0  density.  (cm*3) 

CQBUFF(NS+12 ,NC,NT)  =  H  density.  (cm'3) 


END  of  Standard  MICE/MEGS/MELT  Quantities 
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Optional  Structure  Descriptors 


NQ  =  30  if  ISOURCE  =  1,  2,  or  3  (MICE,  MEGS,  or  MELT). 

=  18  if  ISOURCE  *  4  (SCENARIO) 

If  ISTRIR  =  0,  then  these  structure  quantities  are  omitted. 


CQBUFF(NQ+1,NC,NT) 

CQBUFF(NQ+2,NC,NT) 

CQBUFF(NQ+3,NC,NT) 

CQBUFF(NQ+4,NC,NT) 

CQBUFF(NQ+5,NC,NT) 


=  1st  striation  quantity. 
=  2nd  striation  quantity. 
=  3rd  striation  quantity. 
=  4th  striation  quantity. 
=  5th  striation  quantity. 


See  Table  1. 
See  Table  1. 
See  Table  1. 
See  Table  1. 
See  Table  1. 


END  of  ENVIRONMENT  TUBE  Cell  Quantities 


NOTE:  NTBUFF  is  a  user  computed  implied  DO  loop  limit.  User  must  keep 
track  of  how  many  environment  tubes  have  been  previously  read  from 
the  current  output  block.  If  NTPREV  is  the  number  previously 
read,  and  if  NAVAIL  is  the  number  of  tubes  in  an  output  block 
(JMXBMMXP  if  IFGRID=0,  NATUBES  if  IFGRID=1),  then 

NTBUFF  =  MIN(  NAVAIL-NPREV,  NBUFF/(NQPCEL*JMXA+3)  )  (*) 

This  calculation  must  yield  an  integer  greater  than  zero. 


NOTE:  Proper  interpretation  of  Striation  Quantities  in  Table  1  requires 
background  knowledge  of  methodologies  for  computing  descriptors  of 
plasma  striations.  References  3  and  4  are  recommended. 
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TABLE  1.  Striation  Quantities  Definition. 


For  ISOURCE  -  1  (MICE)  or  2  (MEGS) 


Striation 
Quantity  Number 

|< -  ISTRIP 

1  1  1 

2 

->l 

1 

1 

|  Time  integrated  local  | 

|  growth  rate  j 

Not  Defined 

2 

|  Instantaneous  local  | 

j  growth  rate  (sec'1)  j 

Not  Defined 

3 

|  Not  Defined  | 

1  1 

Not  Defined 

4 

|  Not  Defined  | 

1  1 

Not  Defined 

5 

|  Not  Defined  | 

1  1 

Not  Defined 

For  ISOURCE  *  3  (MELT) 


Striation 
Quantity  Number 

|<  IS 

|  1 

TRIR  >j 

1  2  | 

1 

|  Field-line-averaged 
j  amplitude  (CM) 

|  Field-line-averaged 
|  amplitude  (0—1 ) 

2 

|  a  from  regrids 

1  (cm'3) 

|  <7n  from  regrids 

j  and  convection  (cm'3) 

3 

|  Not  Defined 

1 

1  High n 

i  (cm'3) 

4 

|  Not  Defined 

1 

|  Low  ne 

1  (cm'3) 

5 

|  Not  Defined 

1 

|  High  /  Low 

|  Area  Ratio 
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TABLE  1.  Striation  Quantities  Definition  (continued). 


For  ISOURCE  *  4  (SCENARIO) 


Striation 
Quantity  Number 

(<- 

1 

-  1ST 

1 

RIR  - >| 

1  2  | 

1 

1 

1 

1 

a.  from  convection 
(cm'3) 

[  Reserved  for  algorithm 
j  with  outer  scale  calc. 

2 

1 

1 

Field-mapped  outer 
scale  (cm/rad) 

|  Reserved 

3 

I 

1 

High  ne 
(cm'3) 

|  Reserved 

4 

1 

1 

Low  n_ 

(cm'3) 

|  Reserved 

5 

1 

1 

High  /  Low 

Area  Ratio 

|  Reserved 
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2.7  BETA  FLUX  RECORD  (Optional). 


( ( BETAF ( I , J ) ,  1=1 , JMXB) ,  J=1 , JMXP) 
where: 

BETAF(I.J)  =  Beta  particle  energy  flux  (ergs/sec/cm2) 
through  a  surface  at  100  km  altitude. 

The  surface  in  question  is  contained  entirely 
within  the  (I,J)th  environment  tube  and  is 
perpendicular  to  the  geomagnetic  field.  It  is  the 
center  of  this  surface  which  is  at  100  km  altitude. 


NOTE:  This  record  is  optional.  The  Header  Record  parameter  IBETAR 
indicates  if  this  record  is  present. 
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APPENDIX 


Active  Environment 
Aligned  a  Array: 

Code  Output  Block: 

Environment  Tube: 


GLOSSARY  OF  TERMS 


Tube:  See  "Environment  Tube". 


Boundaries  between  cells  along  the  a  direction 
within  an  environment  tube  (magnetic  flux  tube) 
are  contiguous  with  a  boundaries  in  adjacent 
tubes.  Consequently,  a  cell  boundaries  form 
continuous  grid  surfaces  which  are  everywhere 
perpendicular  to  B.  See  also  "Staggered  a  Array". 


A  cluster  of  records  within  a  file  which  describe 
the  environment  at  a  specified  output  time  in  a 
nuclear  burst  simulation.  Sections  2.1  through 
2.7  describe  the  contents  of  a  code  output  block. 
Different  code  output  blocks  describe  conditions 
at  different  times  or  under  different 
circumstances.  Also  referred  to  as  an  Output 
Block. 


For  output  purposes,  computational  cells  aligned 
with  a  selected  coordinate  direction  have  been 
grouped  in  linear  arrays.  Each  linear  array  can 
be  thought  of  as  a  "tube"  of  cells.  Collectively 
all  of  the  environment  tubes  in  a  code  output 
block  completely  cover  the  three-dimensional 
computational  volume  (subject  to  symmetry 
restrictions  —  see  "Grid  Symmetry"). 

Environment  tubes  from  MEGS,  MELT,  and  SCENARIO 
are  aligned  with  the  ambient  geomagnetic  field 
and,  thus,  constitute  magnetic  flux  tubes.  Each 
cell  in  this  type  of  environment  tube  contains  an 
equal  amount  of  magnetic  flux,  so  the  cross- 
sectional  area  of  cells  varies  inversely  with  B. 

Environment  tubes  from  MICE  are  oriented  along 
radial  lines  from  the  center  of  the  earth. 
Consequently  their  cross-sectional  area  varies  as 
1/R2. 

Historically,  environment  tubes  have  been  known  as 
"sticks",  "tubes",  "flux  tubes",  "plasma  tubes", 
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Grid  Symmetry 


High-#  Plasma 


High/Low  Area 


High  ne: 


or  "columns"  depending  on  context.  Environment 
tubes  may  come  in  two  varieties  —  active  and 
inactive.  Active  tubes  contain  cell  quantities 
which  deviate  from  ambient  due  to  the  action  of 
nuclear  explosions.  Inactive  tubes  contain  only 
ambient  atmospheric  densities,  temperatures,  etc 


No  symmetry  implies  the  complete  three-dimensional 
grid  is  provided  in  the  code  output  block. 

Magnetic  east-west  reflection  symmetry  means 
reflection  symmetry  about  the  vertical  magnetic 
meridian  plane  through  the  center  of  the 
computational  volume  was  imposed  on  the 
simulation.  In  this  case,  only  the  half-space 
used  for  computation  has  been  written  out  in  the 
code  output  block. 

Magnetic  north-south  reflection  symmetry  is 
similar  to  east-west  reflection  symmetry  except 
the  vertical  reflecting  plane  is  oriented 
east-west.  In  MEGS,  MELT,  and  SCENARIO,  the 
reflection  plane  is  the  magnetic  equator.  In 
MICE,  the  reflection  plane  occurs  at  X  =  0. 


Plasma  with  energy  density  well  above  the  magnetic 
energy  density.  Plasma  p  is  the  ratio  of  plasma 
energy  density  to  magnetic  energy  density.  Either 
plasma  thermal  energy  density  or  the  sum  of 
thermal  +  kinetic  energy  densities  is  commonly 
used  in  computing  plasma  p.  Plasma  p  should  not 
be  confused  with  the  p  coordinate  (geomagnetic 
dipole-field-aligned  coordinates  used  in  MELT, 
MEGS,  and  SCENARIO). 


Ratio:  Ratio  of  transverse-to-B  areas  occupied  by  High 

and  Low  ne’s  within  a  computational  cell.  See 
Reference  4  and  "High  ne",  "Low  ne". 


Larger  of  two  electron  densities  assumed  by 
Stagat’s  two-level  closure  relation  in  Striation 
Convection  theory.  See  Reference  4. 


24 


Inactive 
Low  ne: 

MEGS: 

MELT: 


MICE: 


Environment  Tube:  See  "Environment  Tube". 


Smaller  of  two  electron  densities  assumed  by 
Stagat’s  two-level  closure  relation  in  Striation 
Convection  theory.  See  Reference  4. 


A  first-principles,  three-dimensional,  two-fluid 
magnetohydrodynamics  (MHD)  code  for  computing 
evolution  of  high-0  plasma.  MEGS  is  similar  to 
MICE  in  using  detailed  two-fluid  MHD  equations  to 
treat  high-0  plasma;  the  code  is  similar  to  MELT 
in  using  a  global -scale  geomagnetic  dipole- 
field-aligned  coordinate  system.  The  name  MEGS 
derives  from  MHD  Extended  to  Global  Scale. 


A  first-principles,  three-dimensional,  two-fluid 
magnetohydrodynamics  (MHD)  code  for  computing 
evolution  of  low-0  plasma.  MELT’s  use  of  a 
geomagnetic  dipole-field-aligned  coordinate  system 
permits  tracking  of  air  plasma  (generated  by 
high-altitude  nuclear  explosions)  over  global  - 
scale  distances.  Two-fluid  equations  like  those 
in  MICE  are  used  to  represent  physical  processes. 
However,  following  the  assumption  of  low-0  plasma, 
MELT  treats  the  geomagnetic  field  as  ambient. 
Plasma  transport  across  the  geomagnetic  field  is 
via  E  x  B  drift  where  E  is  a  self-consistent, 
spatially  and  temporally  varying  electrostatic 
field.  The  name  MELT  derives  from  Mixed  Eulerian 
iagrangian  Two-fluid. 


A  first-principles,  three-dimensional ,  two-fluid 
magnetohydrodynamics  (MHD)  code  for  computing 
evolution  of  high-0  plasma.  MICE  uses  a 
specialized  geocentric  spherical  coordinate  system 
to  compute  the  dynamics  of  interpenetrating 
neutral  and  ion-electron  fluids.  Collisions  at 
the  atomic  level  transfer  momentum  and  energy 
between  fluids.  Chemical  reactions  convert 
species  of  one  fluid  to  the  other,  as  appropriate. 
Maxwell’s  Equations  are  the  basis  for  computing 
spatially  and  temporally  varying  electromagnetic 
fields.  The  name  MICE  derives  from  MHD  Implicit 
Continuous  fulerian. 
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Electron  density. 


Output  Block: 
Regrid: 

SCENARIO: 


Staggered  a  Array: 


Striation  Convection: 


See  "Code  Output  Block”. 


The  act  of  transferring  information  (i.e.,  cell 
quantities)  carried  on  one  computational  grid  to 
another.  For  example,  one  may  be  motivated  to 
regrid  a  computation  into  a  larger  grid  when  the 
edge  of  an  expanding  fireball  approaches  the 
existing  grid’s  boundaries.  Sometimes  called  a 
"rezone" . 


A  physics'based,  engineering-level  atmospheric 
effects  code  designed  specifically  to  simulate 
multiple  high-altitude  nuclear  detonations. 
Credibility  of  code  outputs  combined  with 
computational  speed  derives  from  use  in  SCENARIO 
of  truncated  forms  of  the  same  three-dimensional, 
two-fluid  magnetohydrodynamics  equations  used  in 
MICE  and  MELT.  SCENARIO  is  available  to  qualified 
users  under  License  from  the  Defense  Nuclear 
Agency. 


Boundaries  between  cells  along  the  a  direction 
within  an  environment  tube  (magnetic  flux  tube) 
are  not  contiguous  or  aligned  with  a  boundaries  in 
adjacent  tubes.  Consequently,  a  cell  boundaries 
do  not  form  continuous  grid  surfaces.  See  also 
"Aligned  a  Array". 


Standard  deviation  of  electron  density;  i.e., 
square  root  of  electron  density  variance.  See 
"Striation  Convection". 


Computational  algorithm  for  estimating  electron 
density  variance  {anz)  produced  by  the  Gradient 
Drift  Instability.  A  specialized  algorithm  is 
required  because  phenomenology  codes  cannot 
compute  fine-scale  plasma  irregularities 
(striations  with  scale  sizes  as  small  as  -  100  m) 
and  simultaneously  span  the  widespread  volume 
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which  is  disturbed  by  evert  one  nuclear  burst 
(i.e.,  hundreds  to  thousands  of  kilometers  in 
extent).  In  the  context  of  plasma  structure, 
striation  convection  is  sometimes  referred  to 
simply  as  "convection".  See  Reference  4. 
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